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Abstract 

We investigate the cosmological gravitino problem in gauge-mediated super- 
symmetry breaking models, where the gravitino becomes in general the lightest 
supersymmetric particle (LSP). In order to avoid the overclosure of the stable 
gravitino, the reheating temperature of inflation Tr should be low enough. Fur- 
thermore, if the gravitino mass is larger than about 100 MeV, the decay of the 
next-to-LSP (NLSP) into the gravitino may modify disastrously the abundances of 
the light elements predicted by the big-bang nucleosynthesis (BBN). We consider 
the case in which the lighter stau is the NLSP and derive cosmological constraints 
from the BBN on the stau NLSP decay. We obtain a lower bound on the mass of 
stau which is more stringent than the current experimental limit m~ > 90 
GeV for the gravitino mass region m 3 / 2 ~ 5 GeV. This lower bound, together with 
the overclosure constraint on the stable gravitino, gives an upper bound on Tr. We 
find that the reheating temperature can be as high as 10 9 -10 10 GeV for m 3 / 2 — 5- 
100 GeV. 



1 Introduction 



Supersymmetry (SUSY) has been considered as an attractive candidate for physics 
beyond the standard model. Many mechanisms have been proposed for SUSY breaking 
and its transmission to our sector. Among them, the gauge-mediated SUSY breaking 
(GMSB) models [|l], |2| are fascinating since they beautifully solve the problem of flavor 
changing processes inherent in the SUSY standard model. Moreover, GMSB models 
are determined by a few parameters, and thus have high predictability. In general, 
they predict that the gravitino is the lightest SUSY particle (LSP) and the next-to- 
lightest SUSY particle (NLSP) is either the lightest neutralino (almost purely bino B) 
or the lighter stau T\ ||. This mass spectrum gives striking features to the low energy 
phenomenology. Here we focus on the cosmological aspects of GMSB models. 

The light stable gravitincQgives rise to a serious cosmological problem, known as the 
"gravitino problem" [f|, || || . The energy density of the gravitinos which are produced in 
the early universe may exceed the present critical density if the gravitino mass is m 3 / 2 ^ 
1 keV [|J. Even if a primordial inflation is assumed, this problem cannot be solved 
completely since gravitinos are reproduced after the inflation ends. At the reheating 
epoch, gravitinos are produced by scatterings in the thermal bath. It has been pointed 
out in Ref. |7j that gravitinos might also be produced non-thermally at the preheating 
epoch. Although this non-thermal mechanism might dominate the thermal production, 
we will not consider it, since it highly depends on inflation models and we would like to 
have a general discussion. 

Stable gravitinos which are produced thermally do not overdose the universe if the 
reheating temperature of the inflation Tr is low enough. The upper bound on Tr is given 
by |, 1 

f 100 GeV to 1 TeV for 1 keV < m 3/2 < 100 keV 

where denotes the bino mass. Since this upper bound on Tr is so severe for a lighter 
gravitino mass regionP], inflation models which predict high reheating temperatures are 
allowed only for a larger 7713/2 region. 

However, when one considers a heavy gravitino (say, m.3/2 ^ 100 MeV) in GMSB, 
there is another cosmological difficulty associated with the NLSP decay. The NLSP 
decays into the LSP gravitino only through gravitational interaction and its lifetime will 

1 We assume here that the i?-parity is exact and hence the LSP gravitino is stable. 

2 Of course, there is no gravitino problem for m 3 / 2 ^ 1 keV [Q. 
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be comparable to the big-bang nucleosynthesis (BBN) era (t ~1-10 2 sec). The decay of 
the NLSP during or after the BBN is dangerous, since the decay products might alter 
the abundances of the light elements, which spoils the success of the BBN || §, [UJ. 
To avoid this difficulty, the lifetime and the energy density of the NLSP are severely 
restricted |, [TL|, 0, M ■ 



In this letter, we investigate cosmological difficulties associated with the stable grav- 
itino and derive the most conservative upper bound on the reheating temperature. Es- 
pecially, we consider in detail the cosmological consequences of the NLSP decay into the 
gravitino around the BBN epoch. Similar analysis had been made in Ref. in which 
the BBN constraint is obtained by considering mainly the effects of hadrons produced 
in the NLSP decay. However, the BBN with high energy hadron injection, examined in 
Refs. H |10| , has not completely been settled yet, since there exist some uncertainties^. 
Here, since we would like to obtain the most conservative constraint on the reheating 
temperature, we will focus on the BBN constraint coming from the photo-dissociation 
of the light elements by the NLSP decay. Furthermore, we consider the case when the 
lighter stau T\ is the NLSP. This is because, compared to the B NLSP, the T\ NLSP has 
larger annihilation cross sections and so its abundance when it decays is smaller, which 
results in a weaker constraint on the reheating temperature. 



2 Estimation of f\ NLSP abundance 

First of all, we discuss cosmological evolution of the f± NLSP and estimate its 
abundance when it decays, since the abundance is crucial for obtaining BBN constraints. 
Let us start from a brief thermal history of T\. By the time when the cosmic temperature 
is comparable to the T\ mass (T ~ only the f\ NLSP among the SUSY particles 

is in thermal equilibrium with the standard model particles. When ri's become non- 
relativistic for T ^ m~ ■ , they pair-annihilate into standard model particles and the 
abundance decreases exponentially. At T ~ Tf (Tj: the freeze-out temperature of rif]), 
Ti decouples from the thermal bath and its abundance freezes out. Finally, when the 
Hubble parameter H becomes comparable to the decay rate, T\ decays into the gravitino. 
The evolution of T\ described above can be traced by solving the Boltzmann equation. 

Before writing down the equation, we must take care of the following two facts. The 
first one is that the masses of charged sleptons f 1 , and t\ (the subscript 1 denotes 

3 For example, there exist uncertainties of the experimental data of the hadron scattering processes, 
and also of the statistical treatment of the errors. 

4 In our case, we find Tf/m~ ~ 1/28-1/25 by numerical calculation. 



2 



the lighter mass eigenstates) are almost degenerate in GMSB models. This is because 
they are mostly right-handed, and hence receive at the messenger scale the same soft 
masses which are determined by the gauge quantum numbers. The mass differences 
between them at the weak scale are induced from the renormalization group effects due 
to the leptonic Yukawa couplings y^ (i = e, /i, r) and also from the left-right mixings 
through Hi. Considering the observed masses of charged leptons, we safely neglect y e 
and and hence the masses for Ji\ and ei, m~ and m~, are the same. On the other 
hand, the tau Yukawa coupling y T gives negative contribution to the stau mass when it 
is evolved from the messenger scale to the weak scale, and the left-right mixing between 
the staus also decreases m~. Both effects lead to the fact that m~ is always smaller than 
m /Ii(ei)- This mass splitting Am = m-^^ — m~ plays a crucial role in calculating the f\ 
abundance. If Am is very small, the relic abundance of T\ is determined by not only its 
annihilation processes, but also by the annihilation with Jl\ and e%, i.e., we should include 
the coannihilation effects [Hfl. We have verified numerically that these coannihilation 
effects become significant for Am ^ Tf. In GMSB models, both Am ^ Tf and Am ^ 
Tj are allowedF], so we respect both cases. 

Secondly, if the bino mass mg is close to m~, we should also take into account the 
coannihilation of f\ with B, as well as with ]i\ and e%. However, this effect increases the 
abundance of f±, which results in more stringent upper bound on Tr (see the following 
discussions). Thus, we forbid the coannihilation of T\ with B and restrict ourselves in 
the parameter region 

m 5> m n+ T /- ( 2 ) 
In the actual calculation below, we take the lower bound on as mg > Lira- to be 
conservative. 

Now we are at the point to present the Boltzmann equation including the coannihila- 
tion effects. This equation describes the evolution of the total number density of charged 
sleptons n — n~ + n~* + n~ x + n~* + n~ + , where Hi stands for the number density 



of one slepton species tF\ [14 



dt 
where 



(171 r i 

-3Hn - (a cS v) In 2 - (n eq ) 2 , (3) 



<^>=£<^}^. (4) 



5 The mass difference Am becomes larger as tan/3 increases, where tan/3 = v u /vd (v u and Vd 
denote the vacuum expectation values of Higgs fields which couple to up-type and down-type quarks, 
respectively) . 

6 We distinguish particles from their anti-particles. 
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Here n eq (nf 1 ) is the equilibrium value of n {jij) and v is the relative velocity of particles 
% and j. The bracket denotes the thermal average and er^- is the total annihilation cross 
section of % + j — » X + X': 



X + X'l 



(5) 



where X and X' represent the possible standard model particles. Note that the Boltz- 
mann equation (|3|) also holds for the case where the coannihilations of n with Ji\ and e\ 
do not occur. 

Next we discuss the thermal-averaged cross sections {a e sv) in Eq. (^). The thermal- 
averaged cross sections {a^v) in Eq. (Q) can be expanded in terms of T/m~ [see Eqs. (|) 
- (§) below]. Since the final abundance is determined by (cr cS v) at T ~ Tf and the 
freeze out temperature is typically Tf ~ m~/25, it is sufficient only to consider the 
leading term, i.e., the s-wave cross sections of the sleptons. Note that the s-wave com- 
ponent of the thermal-averaged cross section (<Tijv) is equal to that of <7jj. The relevant 
(co) annihilation channels are 



(I) 
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Here / denotes the ordinary quarks and leptons. 
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We calculate cross sections for all the processes listed in (I)- (IV). The processes (I) 
and (II) represent the annihilation of each slepton species, while (III) and (IV) describe 
coannihilation processes. In the following we briefly explain their features. First of all, 
let us estimate the cross sections for processes (I). Although they have many final states, 
it turns out that most of them give very small contributions to the total cross section. 
Processes which are induced by the tau Yukawa coupling and/or the left-right mixing of 
sleptons are suppressed due to their smallnessf]. Furthermore, in calculating s-wave cross 
sections, derivative couplings of the sleptons vanish. From these considerations, we find 
that 77 and Z7 are the dominant channels^, and (<7~~*t>) is approximately given byf]: 

(a~~ t v) ~ 2a{nrl - 77) = + O (—) . (6) 

Next, we turn to the processes (II), which have only t-channel and ^-channel diagrams. 
Diagrams induced by the exchange of Higgsino and neutral Wino can be safely neglected 
due to the smallness of y T and left-right mixing, respectively. Then, there is only B 
contribution and its cross section is calculated as 



v 167rac m m~ { T \ 

a~~v) ~ em b + q(J_) (7) 

T1T1 / " cos 4 M^~+^|) 2 WJ 



The annihilation channels listed in (III) can be calculated in the same way as 



a ™ v ) ~ cos* 6 w (ml+m 2 s y + ° ) ' (8) 

for m~ = m~. Cross sections for the processes (IV) are also dominated by the B 
exchange. However, the left-right mixings of sleptons are necessary for these processes, 
and so they give only tiny contributions to the total thermal-averaged cross section 
(a e ffv). Cross sections for other processes in (I)-(IV) are calculated similarly. From 

7 The mixing angles ai of sleptons are in general small in GMSB models (sin a; C 1). However, in 
the extremely large tan/3 region, the stau mixing angle a~ can be large as sina— ~ O(0. 1). Although 

(a~~,v^j /o-(tit* — > 77) might become of 0(10) as opposed to Eq. (||) in this case, we find that the 

final upper bound on Tr only increases by less than a factor of 2. 

8 The ratios of the cross sections for the dominant processes 77, Z7, and ZZ are 1.0:0.6:(0.1-0.2). 
The contributions of the other channels are less than a few percents. 

9 Actually, if m~ < mz/2, the Zj channel is forbidden kinematically. Furthermore, if m~ ~ m h o/2 
(m^o: mass of the lightest Higgs boson h°), the annihilation cross section is enhanced by the pole 
contribution of h°. These effects may change the value of Y in Fig. [lj However, we simply assume 
Eq. @, since the mass region of to— < 90 GeV is excluded from the current experimental limit. 
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m [GeV] 

Figure 1: The final f\ abundance Y. The solid line represents the abundance for 
the case with coannihilation and the dashed line represents that for the case without 
coannihilation. 

the above arguments, we conclude that all the annihilation cross sections are almost 
determined by only the masses of the bino and the sleptons. 

Finally, including all the cross sections described above, we solve the Boltzmann 
equation (|3|) numerically Here we take = 1.1 m~ [see Eq. (Q)]. In this case the 
(co) annihilation cross sections take their maximal values, which gives the smallest f\ 
abundance. In Fig. [1] the result is shown in terms of Y, which is defined by Y = n/s = 
Y^iUi/s with the total entropy density of the universe s. Here are some comments. 
We calculate both of the cases with and without coannihilation effects as mentioned 
before. If the mass difference between jSi(ex) and f\ is large enough (Am^Tj), Jxi and 
Si decay into f\ and disappear from the thermal bath before f\ freezes out. Therefore, 
the coannihilation processes become ineffective. In this case, the abundance becomes 
Y ~ (n~ +n~*)/s for T^zTf, since n~ x — n~ ~ 0. For T^zTj the number density of 
fi and the entropy density decrease at the same rate as the universe expands, and Y 
takes a constant value until f\ decays. In fact, Fig. [I] shows the value of Y just before 
the decay of fx- On the other hand, when Am^Ty, fix an d &i are still in thermal 
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equilibrium at T ~ Tf due to the effect of inverse decays, and thus we should consider 
the coannihilations. Note that Y takes a constant value after sleptons decouple from the 
thermal bath even in this case. This is because Y is an invariant parameter against the 
expansion of the universe and also because n = n~ + n~, + n~ x + n~» + n,~ + n~* does 
not change by the decays of Jli and Z\. Therefore, the final abundance for T\ is given by 
Y for both cases with and without coannihilation effects. 

As shown in Fig. [T], the final abundance with coannihilation is larger than that without 
coannihilation. This can be understood as follows: Imagine when the coannihilation cross 
sections in (III) and (IV) are extremely large. In this case, the final abundance Y would 
be smaller than that without coannihilation. On the other hand, if the coannihilation 
cross sections are zero, Y increase by a factor of 3 because the relevant degrees of freedom 
is now 6, not 2. In our case, coannihilation cross sections in (III) and (IV) are smaller 
than those in (I) and (II), and the final abundance with coannihilation is about twice as 
large as that without coannihilation. 

We present in Fig. [l| the results for the two extreme cases Am ~ and Am ^> Tf, 
for illustration and for comparison. In fact, we find that the final abundance of stau for 
the case ^ Am ^ Tf falls between the two lines in Fig. |l[ 

Before closing this section, we should comment on the case in which the mass differ- 
ence between f\ and /ii(ei) is extremely small. If the mass difference is smaller than the 
tau mass (of course we should include the coannihilation effects in this case), the decay 
channel /xi(ei) — > rxT/i(e) is kinematically forbidden. Then, Ji\ and e\ may dominantly 
decay into gravitinosf^and so we should consider the effect of the Jl\(ei) decay on the 
BBN, as well as t\. However, we simply neglect this possibility, since the final results do 
not change much. 



3 BBN constraint on f\ NLSP decay 

We are now ready to investigate the cosmological consequence of the decay of the 
•fx NLSP into the gravitino at the BBN epoch. The T\ decay rate is estimated, for 

m~ > m 3 / 2 , as 

1 m~ 

T1 " 48vrm2 /2 M2' K ) 

10 The left-right mixing allows Jli to decay through jui — > v^tiT> t . If this is the main decay channel, 
the following discussion does not change, since Jli has already decayed into T\ before the ri decays. 
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Figure 2: Various bounds on the f\ mass. The lower bounds on m~ from the BBN 
photo-dissociation effect are shown by the thick solid and dashed lines for the cases 
with and without the coannihilation effects, respectively. The experimental lower bound 
on m~ is shown by thin solid line. We also present the upper bounds on m~ from 
the constraint Q^J^ p h 2 < 1 by the thick solid and dashed lines for the cases with and 
without the coannihilation effects, respectively. 

with M* = 2.4 x 10 18 GeV, and the lifetime is given by 



n 



6 x 10 4 sec f f ^ \\ (10) 
VI GeV J V100 GeV J V ; 



Therefore, f\ is found to be a long-lived particle and to decay during or even after the 
BBN epoch when m 3 / 2 ^ 10 MeV for m~ = 100 GeV, and hence we should seriously 
consider the effects of its decay on the BBN. 

The energetic r's produced by the f\ decay trigger the electro-magnetic (EM) cascade 
processes and induce high energy photons. These photons may be abundant enough to 
destroy or overproduce various light elements (D, 3 He, 4 He, etc.) synthesized by the 
BBN. In order to keep the success of the BBN, the energy density (per the entropy 
density) of the extra EM particles emitted after the BBN era (i.e., t> 10 4 sec) is severely 
constrained [Q . Here it should be noted that not all the energy of tau contribute to these 
photo-dissociation processes. This is because the tau decays before it causes the EM 
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cascade processes and the produced neutrinos give no effects on the photo-dissociation 
processes 0- By using the result of the recent analysis in Ref. |L2| (see Fig. 16 in Ref. Hl2|), 
we can obtain an upper bound on 

^ = m~Y, (11) 
s 

for a given lifetime of f\. Here denotes the energy density of T\. When one fixes 
the gravitino mass, this bound is translated into the lower bound on m~ , since both the 
abundance and the lifetime of f% are determined by its mass m~ . The obtained result is 
found in Fig. |2|. 

You can see that this BBN photo-dissociation constraint gives a more stringent lower 
bound on m~ for m 3 / 2 ^ 5 GeV, compared to the experimental bound m~ > 90 GeV. 
This result will help us to estimate the upper bound on the reheating temperature in 
the next section. 



4 Gravitino problem and constraint on Tr 

Finally, we discuss the cosmological gravitino problem and obtain an upper bound on 
the reheating temperature of inflation. At the reheating epoch after the inflation ends, 
gravitinos are produced thermally by scatterings with particles in the hot plasma of the 
universe^ 7 ! The relic abundance of the gravitino is given by 0, || 



3/2 ' VI GeV 7 V 100 GeV/ UO 8 GeV, 

where h is the present Hubble parameter in unit of lOOkm/sec/Mpc, and fijw 2 = p|/ 2 / p c 
(P3/2 is the present energy density of the thermally produced gravitinos and p c is the 
critical density of the present universe). It is found from Eq. ([Tj^) that the overclosure 
limit of $7|/ 2 < 1 puts an upper bound on Tr as shown in Eq. (|l|). Notice that this upper 
bound on Tr, if one fixes the B mass, becomes more severe for the lighter gravitino mass 
region, and the highest reheating temperature allowed in GMSB models is achieved for 
a relatively heavy gravitino mass. 

Furthermore, it should be noted that gravitinos are also produced by the f\ NLSP 
decays. Because one gravitino is produced per a f\ decay, we find that 

„nlsp v 
nlsp _ ^3/2 _ ra 3/2 r 

"3/2 = — / ) \ L6 ) 

' Pc Pel S 



11 In fact, the high energy neutrinos produced by the tau decay also causes the EM cascade processes 
by scattering with the background neutrinos. However, this effect is less significant JD|. 

12 As mentioned in Sec. 1, we discard the non-thermal production of gravitinos at the preheating 
epoch in order to obtain the most conservative result. 
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Figure 3: The upper bounds on Tr. The bounds are shown by the thick solid and 
dashed lines for the cases with and without the coannihilation effects, respectively. 



where Py 2 is the present energy density of the gravitinos produced by the T\ decays, 
and so denotes the entropy density of the present universe. In Fig. ^| we also show the 
upper bound on m~ from the constraint Q^jf p h 2 < 1. It is found that the gravitino 
mass is bounded from above as m 3 / 2 ^ 1 TeV when combined with the lower bound on 
m~ from the BBN photo-dissociation effects. 

Therefore, the total relic abundance of gravitinos in GMSB models is now given by 
Q l y 2 = Qy 2 + Q^jf p , as long as the decay of the f\ NLSP takes place within the age 
of the universe. From the lower bound on the f\ mass obtained in the previous section, 
we can estimate the most conservative upper bound on the reheating temperature in 
the following way. Since the abundance Qy 2 depends on the bino mass mg as shown in 
Eq. (|T^), the weakest bound on T R is obtained by the lowest value of mg. On the other 
hand, the lower bound on the mass of the f\ NLSP represented in Fig. ^| is nothing but 
the lower bound on mg. Now we are working in the parameter region of in Eq. (0) 
in order to forbid the coannihilation of T\ with B. Therefore, the lowest value of 
coming from the m~ constraint gives the most conservative upper bound on the reheating 
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temperature^. The result is found in Fig. ^ where we take h = 1 for simplicity^. 

Note that the gravitino mass region of 7723/2 ~ 1 TeV is excluded, since the gravitinos 
produced by the decay of T\ overdose the present universe (^Jf p > 1). On the other 
hand, ^J^ P plays no role in putting an upper bound on Tr for a lighter gravitino mass 
region of m 3 / 2 ^5 500 GeV. 

You can see that the upper bound on Tr is almost proportional to 772,3/2 for 7773/2 ^5 
5 GeV. In this gravitino mass region, the experimental bound on m~ gives the lowest 
value of (see Fig. ||). The upper bound on Tr, thus, is the conventional one Eq.(|l]) 
with ~ 100 GeV. However, when m 3 / 2 ~ 5 GeV-100 GeV, the upper bound on 
Tr takes an almost constant value of 10 9 -10 10 GeV. In this region the bino mass which 
gives the highest value of Tr is determined by the lower bound on m~ from the BBN 
photo-dissociation constraint. Since this lower bound on m~ is more stringent than 
the experimental limit, the upper bound on Tr also becomes more stringent than the 
conventional one. From Fig. 2, the lower bound on the stau mass, i.e., the lower bound 
on the bino mass is almost proportional to (7773/2 for 7773/2 = 5-100 GeV. Thus, it is 
found from Eq.(12) that the upper bound on Tr becomes almost constant. Therefore, in 
GMSB models, the reheating temperature can be taken as high as Tr ~ 10 9 -10 10 GeV 
for 7773/2 ~ 5-100 GeV. 

5 Conclusions and discussions 

In this letter, we have considered the cosmological gravitino problem in GMSB models 
with the f\ NLSP. Especially, we have investigated in detail the cosmological consequence 
of the fi decay soon after the BBN epoch. Since the f\ abundance when it decays is 
crucial for this discussion, we have solved numerically the Boltzmann equations for both 
cases with and without coannihilation effects, and obtained the final T\ abundance. We 
have found that the BBN constraint on the T\ abundance is translated into the lower 
bound on m~, and that the obtained bound is more stringent for 7773/2 ^ 5 GeV than the 
current experimental limit m~ > 90 GeV. This gives some hints for GMSB models. If T\ 
was detected at future collider experiments, the observed f\ mass would enable us to set 
an upper bound on the gravitino mass, and hence on the SUSY breaking scale. It has 
also been found that the gravitino mass region of 7773/2 ^ 1 TeV is excluded (although 
it might be marginal considering the SUSY flavor problem), since the energy density of 
gravitinos produced by the decays of fi's exceeds the present critical density. 

13 In fact, such a value of m~ makes the t\ abundance smallest. (See discussions in Sec. 2.) 

14 The upper bound on Tr becomes more stringent for a smaller h. 
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By using the lower bound on m~ , we have obtained an upper bound on the reheating 
temperature Tr of inflation in order to avoid the overclosure problem of the LSP gravitino 
produced in thermal scatterings and also in the T\ decay. What we have found is that 
the most conservative upper bound on T R in GMSB models is T R < 10 9 -10 10 GeV when 
m 3 / 2 — 5-100 GeV. This upper bound on Tr is weaker than those in the conventional 
hidden sector SUSY breaking models to solve the cosmological problem of unstable 
gravitinos (T R < 10 6 GeV and 10 8 GeV for m 3/2 ~ 100-500 GeV and 500 GeV-1 TeV, 
respectively JT^)- Therefore, it helps us a lot to build SUSY inflation models without 
the cosmological gravitino problem. 

Such a high reheating temperature is also promising for baryogenesis. The relevant 
example is the leptogenesis [16] via decays of heavy Majorana neutrinos, which is very 
attractive from the viewpoint of the observed tiny neutrino masses. The leptogenesis 
requires a high reheating temperature to generate the observed baryon asymmetry of 
the universe^. Thus, the upper bound on the reheating temperature Tr ^ 10 9 -10 10 GeV 
obtained in this letter ensures some scenarios of the leptogenesis to work without the 
gravitino problem. 

The results we have found here are almost independent on models of the GMSB, since 
the abundance of the t\ NLSP is fixed by the f\ mass. However, there are some loop- 
holes in which these constraints can be avoided. If one introduces the .R-parity breaking, 
the NLSP can decay before the BBN (i.e., 1 sec) and evade the BBN constraints, 
and also the LSP gravitino can decay within the age of the universe. Furthermore, if 
one assumes the late-time entropy production such as the thermal inflation [19| in the 
thermal history of the universe, the abundances of the NLSP and also the LSP gravitino 
are diluted away so that we are free from these constraints. 

In the present analysis, we have considered only the photo-dissociation effects of the 
decay of f\ on the BBN to make a conservative analysis. As pointed out in Ref. [13| the 
high energy hadrons produced by the t\ decay might also be dangerous. Here we would 
like to briefly comment on their effects. The high energy hadrons, if they are produced 
at the time t ~ 1-10 4 sec, delay the freeze-out of the p-n conversion and raise the number 
ratio of n to p, which leads to the overproduction of D and 4 He JTIJ. Furthermore, if T\ 
decays at the time t ^ 10 4 sec, the produced hadrons destroy the light elements synthe- 
sized by the BBN and modify their abundances (e.g., increases the 7 Li abundance) M. 



15 In the leptogenesis scenarios where heavy Majorana neutrinos are thermally produced, the reheating 
temperature of Tr ^ 10 10 GeV is required to induce the sufficient baryon asymmetry JT^ . On the other 
hand, for the case when heavy Majorana neutrinos are non-thermally produced in the inflaton decays, 
the reheating temperature of Tr^ 10 6 GeV is required [fill. 
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These effects give the upper bound on the f\ abundance. However, there are some uncer- 
tainties in the BBN with high energy hadron injection. The relevant hadron scattering 
cross sections, especially those of Li, have not been observed experimentally in detail. 
Furthermore, the statistical treatment for the estimation of the errors has not completely 
settled yet. Therefore, in the present analysis we do not include the BBN constraints 
from the hadron injection. 
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